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Abstract 

In this paper, we consider nonlinear diffusion processes driven by space-time white 
noises, which have an interpretation in terms of partial differential equations. For 
a specific choice of coefficients, they correspond to the Landau equation arising in 
kinetic theory. A particular feature is that the diffusion matrix of this process is a 
linear function the law of the process, and not a quadratic one, as in the McKean- 
Vlasov model. The main goal of the paper is to construct an easily simulable diffusive 
interacting particle system, converging towards this nonlinear process and to obtain 
an explicit pathwise rate. This requires to find a significant coupling between finitely 
many Brownian motions and the infinite dimensional white noise process. The key idea 
will be to construct the right Brownian motions by pushing forward the white noise 
processes, through the Brenier map realizing the optimal transport between the law 
of the nonlinear process, and the empirical measure of independent copies of it. A 
striking problem then is to establish the joint measurability of this optimal transport 
map with respect to the space variable and the parameters (time and randomness) 
making the marginals vary. We shall prove a general measurability result for the mass 
transportation problem in terms of the support of the transfert plans, in the sense of 
set- valued mappings. This will allow us to construct the coupling and to obtain explicit 
convergence rates. 

Key words and phrases: Landau type interacting particle systems, nonlinear white noise 
driven SDE, pathwise coupling, measurability of optimal transport, predictable transport 
process. 
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1 Introduction and main statements 

Consider the nonlinear diffusion processes in M"^ of the following type: 

Xt = X^+ f I a{Xs - y)Wp{dy, ds) + f I b{X, - y)Ps{dy)ds (1) 
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where Pt is the law of Xt, and Wp is a valued space-time white noise on [0,T] x M.'^ 
with independent coordinates, each of which having covariance measure Pt{dy) (8) dt. 

The nonlinear process ([1]) was introduced by Funaki [3j , who obtained existence and unique- 
ness results for Lipschitz coefficients a : M.'^ ^ ■^d(g)d b : M.'^ ^ M.'^ , see also Guerin [7] 
for a different approach. It has an important interpretation in terms of partial differential 
equations issued from kinetic theory. More precisely, for a specific choice of coefficients a 
and b, the laws {Pt)t are a weak solution of the spatially homogeneous Landau (also called 
Fokker-Planck-Landau) equations for Maxwell potential: 
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df df 

f{t,v.) — {t,v)-f{t,v)g^^{t,v. 



dv^ 



(2) 

with aijiy) := {aa*)ij{v) = \v\'^5ij — ViVj and bi{v) = V • aj.(f). The equations ([2]) model 
collisions of particles in a plasma and can be obtained as limit of the Boltzmann equations 
when collisions become grazing, see Funaki [4J, Goudon [5], Villani [T?j ^18j and Guerin- 
Meleard [8]. 

In this work, we shall prove the convergence in law of an easily simulable mean field inter- 
acting particle system towards the nonlinear process ([T|) at an explicit pathwise rate. This 
problem is of great interest in order to construct a tractable simulation algorithm for the 
law Pt and thus, in particular, for solutions / of equation ([2]). To our knowledge, there 
is no result on convergence rates of the deterministic numerical methods used at present 
for the Landau equation, which are reviewed in [2]. The interest of our approach is that 
it is based on the diffusive nature of the equation, and that it addresses a large class of 
nonlinear processes. The fact that we want to deal with simulable systems will necessitate 
a coupling between finite dimensional and infinite dimensional stochastic processes. We 
shall introduce a coupling argument based on new results on measurability of the optimal 
mass transportation problem. 

We consider a particle system which is naturally related to the nonlinear process. Indeed, 
notice that the diffusion matrix associated with ([l]) is defined on M'^ by 

a{x,Pt):= [ a{x-y)a*{x-y)Pt{dy) = [{aa*)*Pt]{x). (3) 

Thus, if in order to approximate the white noise driven stochastic differential equation ([T]), 
we heuristically replace Pt in ([3]) by an empirical measure of n E N* particles in M'^, we are 
led to consider the following system driven by independent Brownian motions (i?*^): 

Xr = X^, + — / ^a{Xi'--X^^-)dBf + - / ^b{Xl'--X^'-)ds,i = l,...,n. (4) 

To be more precise, if /i" = ^ Y17=i ^x^'" empirical measure of the system, the 

mappings 

1 r* " 

fit, u;,x)^^ V /(s, a;, X^^^)dBl\ i=l,...,n, (5) 
-^0 k=l 

define (for suitably measurable functions /) orthogonal martingale measures in the sense 
of Walsh [20], with covariance measure dt. 
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By adapting techniques of Meleard-Roelly based on martingale problems, one can show 
propagation of chaos for system with as limit the process ([1]). This says in particular 
that the covariance measure of ([5]) converges in law to Pt dt when n goes to infinity. But 
in turn, the arguments of [11] do not give any information about speed of convergence. 
To estimate the distance between the law of the particles and the law of the nonlinear 
process, we need to construct a significant coupling between finitely many Brownian motions 
and the white noises processes. This problem is much more subtle than in the McKean- 
Vlasov model (cf. Sznitman [16] or Meleard [E]), where each particle is coupled with a 
limiting process through a single Brownian motion that drives them both. The well known 
convergence rate in that model is consequence of the standard L^-law of large numbers 

in and of the fact that the diffusion and drift coefficients of the nonlinear process depend 
linearly on the limiting law through expectations with respect to it. In the present Landau 
model, we have to deal with the space-time random fields ([5]), which have fiuctuations of 
constant order in n. This is also refiected in the fact that it is the squared diffusion matrix 
of that depends linearly on Pt (see ([3])). It is hence not clear where a convergence rate 
can be deduced from. 

Let ,i = 1, . . . ,n he n independent copies of the nonlinear process in some probability 
space, and their empirical measure at time t (observe that it samples Pt). We shall 
construct particles ([!]) on the same probability space, in such way that they will converge 
pathwise in on finite time intervals, at the same rate at which the Wasserstein distance 
W2 between Pt and goes to 0. Let us state our main result on the process ([1]): 

Theorem 1.1. Let n G N and assume usual Lipschitz hypothesis on a and b, and that the 
law Pq of Xq has finite second order moment. Assume moreover that Pt has a density with 
respect to Lebesgue measure for each t > 0. 

Then, in the same probability space as {X^, . . . , X^) there exist independent standard Brow- 
nian motions (-B**^)i<i^fc<„ such that the particle system (^*'"')"=i defined in ^ satisfies 

e( sup \Xr-XlA <CeMC'T) E{W^{u:,P,))ds 
\te[o,T] J Jo 

for constants C, C that do not depend on n. 

Thanks to available convergence results for empirical measures of i.i.d samples (see e.g. 

|14j). Theorem 1 1 . 1 1 will allow us to obtain, under some additional moment assumptions on 

-2 

Pq, the speed of convergence for the pathwise law of the system (see Corollarv 16. 2p . 
We remark that the absolute continuity condition of Theorem 11.11 can be obtained under 
non-degeneracy of the matrix aa* by using for instance Malliavin calculus [13j; it is also 
true for the specific coefficients of the Landau equation ([2]) despite their degeneracy, and 
for some generalizations (see Guerin [6j). 

The proof of Theorem 11.11 relies on new results on the optimal mass transportation problem. 
For general background on the theory of mass transportation, we refer to Villani |19j . 
Recall that if /i and u are probability measures in with finite second moment, the first 
of them having a density, then the optimal mass transportation problem with quadratic 
cost between ^ and v has a unique solution, which is a probability measure on M^'^ of the 
form Tr{dx,dy) = ii{dx)5j'(x)[dy) . The so-called Brenier or optimal transport map T{x) is 
(// a.s. equal to) the gradient of some convex function in W^, and pushes forward /i to v. 
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Let now Wp be the white noise process driving the z-th nonhnear process X*. The key 
idea in Theorem 11.11 wiU be to construct Brownian motions (-B*'^)fc=i...ra in an "optimal" 
pathwise way from Wp. Heuristicahy, this wih consist in pushing forward the martingale 
measure Wp through the Brenier maps T*''^'"'(x) realizing the optimal transport between 
Pt and v'^ioj) (this is the reason for the absolute continuity assumption on Pt). But to give 
such a construction a rigorous sense, we must make sure that we can compute stochastic 
integrals of T*''^'"(x) with respect to Wp{dx.,dt). Prom the basic definition of stochastic 
integration with respect to space-time white noise (cf. |20j). this requires the existence of a 
measurable version of (t,u;,x) i— > r*''^'"(a;) being moreover predictable in (f,u;). A striking 
problem then is that no available result in the mass transportation theory can provide 
any information about joint measurability properties of the optimal transport map, with 
respect to the space variable and some parameter making the marginals vary. Nevertheless, 
we will show that a suitable "predictable transportation process" exists: 

Theorem 1.2. There exists a measurable process {t,uj,x) i— > T'^(t,uj,x) that is predictable 
in {t,u}) with respect to the filtration associated to {Wp, . . . ,Wp) and (Xq,...,Xq), and 
such that for dt ¥(duj) almost every {t,uj), 

T"-{t,LO,x) = r*''^'"(x) Pt{dx)-almost surely. 

This statement is consequence of a general abstract result about "measurability" of the mass 
transportation problem. To be more explicit, recall that the optimality of a transfert plan vr 
is determined by its support (it is equivalent to the support being cyclically monotone, see 
McCann |lOj or Villani [H]). On the other hand, without assumptions (besides moments) 
on the marginals /i and v, the solution vr of the mass transportation problem may not 
be unique. A basic question then is how to formulate, in a general setting, the adequate 
property of "measurability" of the solution(s) vr with respect to the data (//, v). As we shall 
see, the natural formulation requires to introduce notions and techniques from set-valued 
analysis. Then, we shall prove the following 

Theorem 1.3. Let 7^2 (I^"^) be the space of Borel probability measures in with finite 
second order moment, endowed with the Wasserstein distance and its Borel a— field. Denote 
by n*(/i,i/) the set of solutions of the mass transportation problem with quadratic cost 
associated with (/i,z^) G {V2(^'^))'^ ■ The function assigning to {n,v) the set o/M^°'.' 

supp(tt), 

is measurable in the sense of set-valued mappings. 

In particular, this ensures that if and vx vary in a measurable way with respect to some 
parameter A, so that in each of the associated optimal transportation problems uniqueness 
holds, then the support of the solution vr^ also "varies" in a measurable way. This will be 
the key to our results. 

The rest of this work is organized as follows. In Section 2 we review the Wasserstein 
distance and the mass transportation problem with quadratic cost in (in particular the 
characterization of its minimizers) . In Section 3 we prove Theorem 11.31 and a consequence 
needed to prove Theorem 11.11 In Section 4, we state some properties about process ([1]) 
and we heuristically describe our coupling between space-time white noises and Brownian 
motions. In Section 5 we construct the "predictable transportation process" of Theorem 
11.21 needed to rigorously define the coupling. Section 6 is devoted to complete the proof of 
Theorem 1 1 . 1 1 and to obtain explicit convergence rates. 
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2 The mass transportation problem with quadratic cost in 
and the Wasserstein distance 

We denote the space of Borel probability measures in R*^ by 7^(M'^), and by 7^2 (IK'') the 
subspace of probabihty measures having finite second order moment. 

Given vr G V2(^'^'^), we respectively denote by tti and tt2 its first and second marginals on 
M'^. On the other hand, for any two probability measures ^, G V2{^'^) and vr G 
we write 

if TTi = /u and 7r2 = i'. Such vr is refereed to as a "transfert plan" between fi and ly. 
Definition 2.1. The Wassertein distance W2 on V2{R'^) is defined by 



z^) := inf / \x — y\'^7r{dx,dy). 



Then, ^2^ j^'^'i) is a Polish space, see e.g. Rachev and Riischendorf [14J. The topology 
is stronger that the usual weak topology. More precisely, one has the following result (see 
for instance Villani, |19] Theorem 7.12) 

Theorem 2.2. Lei G 'P(M'^). The following are then equivalent: 
i) W2{fJ-^,^J') — > when n — > 00. 
a) /i" converges weakly to ji and 



\x\ fj,"'{dx) I \x\ fj,{dx 

Hi) We have 

ip{x)f/^{dx) — > / ip{x)fi{dx) 



for all continuous function c/? : M"^ ^ M such that \^{x)\ < C(l + |xp) for some C G 
We shall denote by L the mapping L : 7^2(1^2'^) ^ M defined by 

L{it) = / \x — y\'^TT{dx, dy). 



Remark 2.3. It is not hard to check that L is lower semi continuous (l.s.c) for the weak 
topology. Moreover, L is continuous for the Wasserstein topology in V2(^'^'^) by part Hi) of 
Theorem \2.2l 



Fix now /i, 1/ G 7^2 (M''), and denote by n*(/i, u) the subset of 7^2 (K^"*) of minimizers of the 
Monge-Kantorovich transportation problem with quadratic cost for the pair of marginals 
{n, u) . That is, 

n*(/i,i/) := argmin^^^L{'K). 

It is well known that 11* (/i, v) is non-empty. Indeed, it is not hard to see that for the weak 
topology, {tt G 7^2 (IK^*^) : <v} is a compact set, and the lower semi-continuity of L implies 
the existence of minimizers (see e.g. |19] Chapter 1 for details). 

We shall next recall the characterization of minimizers of the transportation problem with 
quadratic cost. We need the notion of sub-differential of a convex function: 



5 



Definition 2.4. Let : A d ^] - 00,00] he a proper (i.e. if ^ +00J lower semi- 
continuous (l.s.c) convex function. The suh- differential of ip at x is 

dip{x) = {y eR"^ : (f{z) > ip{x) + {y,z - x),yz e M'^}. 

Elements of d(p{x) are called sub-gradients of ip at point x. The graph of dip is 

Gr{dip) = {{x,y) G M^"^ : y G d^{x)] 

and it is a closed set. 

Recall that (p is differentiable at x if and only if d(p{x) is a singleton (in which case d(p{x) = 
{V(p{x)}). Also, the set {x G M*^ : 93 is differentiable at x} is borelian, see e.g. McCann 

m- 

We next summarize results in pioneer works in this domain, Knott-Smith [9], Brenier [1] and 
McCann [lOj . Rachev and Riischendorf [14]. See also Villani [19] for a complete discussion 
on these questions, proofs and background. 

Theorem 2.5. Let fj.,u e V{R'^) and IT <u be a transfert plan. We have 

a) TT G n*(M) 1^) if o-nd only if there exists a proper l.s.c. convex function 99 such that 

supp{-k) C Gr{dip) 

or, equivalently 

7r{{{x,y)eM.^ :yedip{x)}) = l. 

b) Assume that jj, does not charge sets of Hausdorff dimension less or equal than d — 1 
and that vr G 11* Then, 

i) the set {x G M'^ : 99 is not differentiable at x} has null ^-measure, 
ii) We have 

TT{dx, dy) = n{dx) ® 8^^(x){dy)- 

ii) If T is a measurable mapping such that TT{dx,dy) = fi{dx) Sx(x){dy), then 
T{x) = Vip{x) , fj,{dx) — a.s.. 

Hi) vr G n*(/i,z/) is unique. 

This result will be useful later in the particular case when the measure /i is absolutely 
continuous with respect to Lebesgue measure. 

3 Measurability of the mass transportation problem 

We now introduce the basic notions on "multi-applications" or "set-valued mappings" that 
we need to prove Theorem 11.31 For general background, we refer the reader to Appendix 
A in Rockafellar and Wets |15j . 

Definition 3.1. Let X,Y be two sets. 

i) A function S on X taking values in the set of subsets of Y is called a set-valued 
mapping or multi- application. We write S : X . 
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a) For any A<zY , the inverse image of A through S is the set 

S-\a) ■.= {xeX : S{x) n A / 0}. 

Hi) If{X,A) is a measurable space and {Y,Q) a topological space, we say that S : X 
is measurable if for all 9 G Q, 

S~He) G A. 

(Of course, if S{x) = {s{x)} is singleton for allx, measurability of S is equivalent to 
that of s. ) 

Consider 7^2 (I^'^) endowed with the Wasserstein distance and the Borel cr— field. We define 
a set- valued mapping 

by 

^'(/u,z/) := {{x,y) : 3-k G 11* s.t. {x,y) £ supp{tt)}. 

Our goal is to prove that ^ is measurable. We shall need some further notions on set-valued 
mappings. 

Definition 3.2. Let X be a set, and (Y, H) and {Z, 0) be topological spaces. 

i) A set-valued mapping S : X ^Y is closed-valued if for all x £ X, S{x) is a closed 
set of{Y,E). 

ii) A set-valued mapping U :Y ^ Z is inner semicontinous (i.s.c) if for all 6 € Q, 

The following results can be found in Appendix A of [15], in the case of set-valued mappings 
in M'^. For completeness we provide proofs in a more general context. 

Lemma 3.3. Let {X,A) be a measurable space and {Y,E) a topological space. 

i) S : X ^ Y is measurable if and only if the closed-valued mapping x ^ Cl{S{x)) is 
measurable, where Cl{S{x)) is the topological closure of the set S{x). 

ii) Assume that {Y,d) is a separable metric space and that S : X ^Y is closed-valued. 
Then, S is measurable if and only if for all closed set F of Y , 

S-\F) G A. 

Hi) Let {Y, E) and (Z, 0) be topological spaces, S : X ^Y be measurable and U :Y ^ Z 
be i.s.c. Then, the multi- application U o S : X ^ Z, defined by 

UoS{x):= U U{y) 

yeS{x) 

is measurable. 
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Proof i) For any open set 9 e E, S{x) n 6^ 7^ if and only if Cl{S{x)) n 6* / 0. 

ii) "Only if" part: since y is a metric space, we use that every closed set F is the intersection 

of some countable collection of open sets {9n)- Therefore, 

{xeX : S{x) nF/0}= f]{xeX : S{x) n 0n / 0} G A 

neN 

"If" part: {Y, d) being separable, we can express every open set 9 as the union of some 
countable collection {Bn) of closed balls. We then have that 

{x^X: S{x) n0/0}= U{xGX: S{x) n S„ / 0} G ^. 

neN 

in) Straightforward: 

{U o S)-\9) = {xeX: {Uy^s{x)U{y)) n9^$} = {xeX : 3y e S{x) s.t. U{y) n M 0} 
= {xeX:S{x)r\{U-^{9))^$]. 

The function U being i.s.c, U^^{9) belongs to H, which allows us to conclude. 

□ 

Now we can proceed to the 
Proof of Theorem 11.31 

We observe first that ^'(/i, zv) = U o zv), where S and U are the set valued mappings 
respectively defined by 

and U : 7^2(^2'^) ^ R'^ by 

U{tt) := supp{'k) 

We will therefore split the proof in several parts: 

a) S is a dosed valued mapping 

First notice that tt 1— > Tr^ is continuous for the Wasserstein topology. Indeed, VF2(7r"', vr) — > 
implies that tt" converges weakly to vr, and then vr" converges weakly to tTj for i = 1,2. More- 
over, we have f^a \x\'^iriidx) = f^2d \x\'^Tr^{dx,dy) J^2d \x\'^TT{dx,dy) = J^a \x\'^Tri{dx) by 
Theorem 12.21 and then the asserted continuity follows. 
Consequently, vr 1-^ W2{iti,tt2) too is continuous. Therefore, 

n* (//,!/) = : vr n {tt : L{tt) - ^2(^1, ^2) = 0} 

is the intersection of two closed sets 7^2 (IR^'^)- 

b) Inverse images through S of closed sets are closed sets 

Let F C V2iR'^) be a closed set and G S'^F), n G N , be a sequence converging 

to (/i, u) in {V2(R'^))'^ ■ Then, ^ ^ and u"^ ^ u weakly, and (^") and {u"') are tight. 
But since (/i",!^"") G S~^{F) for each n, there exists iTn s.t. vr" <(^n, and then (7r„) too is 
tight (by considering products of compact sets) . 

Let (tt"*) be a weakly convergent subsequence with limit tt. Then, clearly vr <u. We will 
prove that L{-k) = W2{fJ',i^) and that tt £ F, which will mean that (/i, z^) G S~^{F) and 
finish the proof. 
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We have 

/ (|a;|2 + |y|2)vr"'=(dx,d2/) = / |x|V=(d3;)+/ |ypi/"'=(dy) ^ 

/ \x\'^n{dx)+ \y\^u{dy)=l (Ixp + |y|^) 7r(dx,(iy), 

which imphes that W2{tt'^,tt) — > and vr G F. Finally, by the continuity of vr i— > -L(7r) — 
W2{iTi,'7T2) we get that 

= L(7r"'=) - W^sKNvr^'^) = Htt) - W2{fi,u). 

c) The mapping U is i.s.c. 

Let be an open set of M?'^. We must check that 

{vr G V2{M.'^'^) : suppCvr) n ^ / 0} = {vr G ^2(^2°') : 7r(6') > 0} 
is open, or equivalently, that 

{tt G T'2(K2'^) : vr(0) = 0} 

is closed in V2{M.'^'^). Assume that vr.vr" G P2(K^"'), with vr" such that 7r"(6') = for ah 
n G N, and moreover that VF2(7r",vr) ^ 0. Then vr" converges weakly to vr, and so by the 
Portemanteau theorem, we have 

O = liminfvr"(0) > vr(e). 
n 

d) Conclusion 

By parts a) and b) and Lemma 13.31 ii) we get that 5* is measurable. By c) and Lemma 13.31 
Hi) U o 5 is measurable and the proof is finished. 

□ 

The following corollary will be useful in the specific setting needed to prove Theorem ll.il 

Corollary 3.4. Let {E,Tj) be a measurable space, and X G E {l^x,i^x) £ (P2(I^'^))^ and 
^ : E ^M."^ be measurable functions. Then, the set 

{(A,x):(x,e(A))GC/W(^A,^A)} 

belongs to T. (g> B{W^) 

Proof By Lemma 13.31 i) and Theorem 11.31 we get that Cl{^) is measurable. Moreover, it 
is not hard to check that the mapping 

is measurable and closed-valued. Then, we just have to notice that 

(x, e(A)) G Cmifix, i^x) if and only if [Cmifix, ^x) - (x, ^(A))] n C / 
for the closed set C = {0}. 

□ 
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4 A coupling between space-time white noise and Brownian 
motions via optimal transport 



In all the sequel, we refer the reader to Walsh [20] for background on space-time white noise 
processes and stochastic integration with respect to martingale measures. 
Assume that a : R'^ ^ M*'^ and b : ^ are Lipschitz continuous and with linear 
growth. Then, by results of [3] or [7j we can construct in some probability space (r2,J^, P) 
a sequence of independent copies of the nonlinear processes, 

XI = X^o+ f [ a{Xl - y)Wh{dy, ds) + f I h{Xl - y)Ps{dy)ds, (6) 

Jo JRd Jo Jr^ 

where the Wp are independent space-time R'^-valued white noises defined on [0, oo) x M"^. 
Each of the d (independent) coordinates of Wp has covariance measure Pt{dy) dt, where 
Pt is the law of Xt. The initial conditions {Xq, . . . , Xq, . . . ) are independent and identically 
distributed with law Pq, and independent of the white noises. The pathwise law of X* is 
denoted by P, and it is uniquely determined. 

Denote by J-J^ the complete right continuous c-field generated by 

{iW}.{[0, s] X A^), W^i[0, s] X ^")) ■.0<s<t,A' £ BiR'^)} 
and {Xq, . . . , Xq). We also denote by 

the predictable field generated by continuous (.F")-adapted processes. 

In what follows, we fix a finite time horizon T > 0. Under usual Lipschitz assumptions on 
the coefficients, there is propagation of the moments of the law Pq, as proved in Guerin [7]. 

Lemma 4.1. If E[\Xq\^) < oo for some k >2, then 

E sup \Xt\^ I < oo. 
V*e[o,T] ) 

The continuity of X and the previous uniform bound imply that t i— > jj^^j \x\^Pt{dx) is 
continuous. 

Throughout the sequel, the assumptions of Theorem 11.11 on Pq and Pt are enforced, in 
particular, the condition £'(supjg[o,T] l^tP) < ^ '^ill hold by the previous lemma. 

We shall now present the main idea of the coupling we introduce to prove Theorem II. li 
Basically, this consists in constructing for each n, Brownian motions in a pathwise way, 
from the realizations of the n white noises {Wp, . . . , Wp). The key for that will be to 
use the optimal transport maps between the marginal Pt of the nonlinear process and the 
empirical measures of samples of that law. More precisely, write 

1 " 

1=1 
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and notice that for each uj e n, < t < T) is an element of C{[0,T],V2{^'^))- Thus, 

for each i G [0, T], n S N and lo, and we can consider the optimal coupling problem with 
quadratic cost between I't^uj) and Pt, 

inf < / |x — yp7r(dx, dy) > . 

By the assumption on Pt and Theorem 12.51 the following properties hold for each fixed 
pair {t,io) e]0,T] x Q: 

Lemma 4.2. a) There exists a unique tt*''^''^, such that 

WiiPt, u^u^)) = [ \x- y| 

b) There is a Pt{dx) — a.e. unique measurable function T*'"^'" : M'^ — > M'^ such that 

TT^'^''^{dx^dy) = 6rpt,u,,n(^^-^{dy)Pt{dx) . 
In particular, under Pt{dx) the law of T*'^'^{x) is 

c) We have 

WiiPt, u^iiv)) = [ \x- T''^'^[x)\^Pt{dx). 

We would like to construct independent Brownian motions by "transporting" the n 
independent white noises (W^, . . . , VFp) through the transport mappings T*'^'"'(a;). As 
pointed out in the introduction, to do so we must at least be able to define stochastic 
integrals of functions of the form {t,uj,x) ^ /(T*''^'"(x)), with respect to the white noise 
processes. The existence of a version T"(t, u, x) of T*''^'"'(x) having good enough properties, 
will be established in next section, when we shall prove Theorem 11.21 

Before doing so, we observe that if Theorem 11.21 holds . then the following processes Bf' = 
B^^'^ will be well defined from Q. 

Proposition 4.3. For each n G N*, define 

-Sf '"(w) ■■=Vnj [ i{T^{s,c,,x)=xHi^)}^p(dx, ds), i,k = l...n (7) 
Jo JR'' 

Then, (-B*'^'"')i<j,A;<n are independent standard Brownian motions in W^. 

These are the right Brownian motions we need to construct dH). The proof of Proposition 
14.31 will given in Section 6. 

5 Construction of the predictable "transport process" 

Our goal now in this section is to show that for each n G N*, there exists a process 
{t, oj, x) I— > T^{t, u), x) defined P(da;) ® dt® Pj ((ix)-almost everywhere, which is measurable 
with respect to Vred^ ® B(R'^), and such that: 

for dt 0F{du)) almost every 
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r"(t, cj, x) = r*''^'"(x) Pt(dx)-almost surely . 

Since {X^, . . . are independent copies of the nonlinear process and Pt = law{Xl) has 
a density with respect to Lebesgue measure, for each t G [0, T] we have that 

n^i^j:Xl = Xl)=Q. 

Notice also that for fixed (i, j) with i ^ j the following set 

{{t,u):Xl{u;)=Xl{u)} 

belongs to VredP" since (i, cu) i— > \Xl{uj) — Xl {uj)\ is adapted and continuous in t. 
By Fubini's theorem we then see that 

lr^«=^ji(i,w)P(dw) (it = 

[0,T]xn ^ ^ 

Remark 5.1. Consequently, there is a predictable set of [0,T] x 

of full ¥(duj) <Ei dt-measure and such that 

for all {t, iu) £ Qlj,, XI{lv) / Xi{uj) for all i, j G {1, . . . , n}. 
Let us denote by (Vred^)' the a— field Vred^ restricted to il^. 

Recall that for each the set of solutions 11* ^P^, ^ °f ^he optimal trans- 

port problem between Pt and ^ X]r=i ^xi^^) ^ singleton that we have denoted by tt*'"^'". 
Let us define now the sets 

:= {(t,u;,x) G Vl'^ x M'^ : G S'upp(^*''^'")} , i = l,...,n. 

The sets ^4*'" are predictable, as proved in the following lemma. 
Lemma 5.2. We have G (Pred")' ^(M"^). 

Proof Observe that the deterministic process {t,oj) ^ Pt e V2{^'^) is Pred"-measurable. 
Indeed, if (/n)neN\{o,i} is a countable dense subset of the space of continuous functions in 
with compact support, and fo{x) = l,fi(x) = then the topology of P2(I^'^) is generated 
by the real mappings m 1-^ / fn{x)m{dx). It is therefore enough that (t, w) J fn{x)Pt{dx) 
be Pred^-measurable, which is clear since t ^ Pt is continuous. 
Next we will apply Corollary 13.41 to the measurable space 

(i?,S) = (l^^,(Pred")'), 
A = and the (7^red"')'-measurable functions given by 



{t,u) 



Pt, ^ 6^, (w) j G 7^2(^2'^) and {t, Lo) ^ e{t, uj) = Xfiuj) G 
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For each {t,uj) E 0^, with ^ denoting the multi-apphcation defined in Theorem 11.31 we 
simply have in the current setting that 



Corollary 13.41 impHes the result. 

□ 

Recall from basic measure theory that if Ei and E2 are measurable spaces and A C EiX E2 
is an element of their product cr— field, then, for each Ai £ Ei, the fiber of A at Ai is the 
set 

[A]x, := {XeE2:{\i,X)eA}, 

and it is always measurable in E2. 
We can now proceed to the 

Proof of Theorem II. 2t 

We split the proof in several parts. 

a) The sets A*'", i = 1 . . .n form a partition of Jl^ x up to F{duj) (E>dt(E> Pt{dx)-null sets. 
For i ^ j write 

:= {{t,uj,x) G (fi^ X M'^) : {x,Xl{uj)) G S'upp(7r*''^''^) and {x,Xi[uj)) G supp(7r*''^''^)} 



and denote by := {x e : it,uj,x) G A^^-"} G ^(M'^) the fiber of A^J> at 

(t,tj) G 0,'rp. Then, we have 

mA'^M) = e : ix,Xliu;)), {x,Xl{u,)) G 5npp(7r*'-'")}) 

< Pt{{x G : Xi{uo),Xi{uo) G 

where f^*'"^'" is a proper l.s.c. convex function given by Theorem 12.51 a). But since (t, cj) G 
Jl^, we have XKu) 7^ and so 

Xl{u)),xi{Lo) G =^ if is not differentiable in x. 

We obtain by Theorem [23] that ( (j^^) ) = 0, and then 

e \ I l^i,„n^j,n(t,u;,x)Pt(dx)dt =0. 

yj[0,T]xnxR'* j 

On the other hand, since r*''^'"(a;) G {^^^(a;), . . . ,X"(a;)} Pt{dx) a.s., we have for all (f, w) 
that 



Pi 



\i=l 



= G M'^ : for alH = 1, . . . , n, (x, supp(7r*'^'")}) 

< G M'^ : (x,r*'^'"(x)) supp(^*''^'")}) 

= 7r*'"'"(sMpp(7r*'^'")^) 
= 
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Defining the set 



we deduce that Pt \ \Q.^](t^^-^j = for all {t,uj) £ 17^. Therefore, 



and so fir has full P((iw) iE> dt<Si Pt{dx)— measure. This proves assertion a). 
We can now define a Vred^ Cg) ;S(M'^)— measurable function by 

n 

r"(t,u;,x) := j;i^„„nn^(t,^,x)XKa;). (8) 
1=1 

6j For P(dLij) (8) dt almost every (t,uj), T^{t,uj,x) = T*''^'"(x) /lo/c/s Pt{dx) almost surely. 
By Theorem 12.51 this is equivalent to prove that 

7r*''^'"((ix, dy) = Pt{dx) (g) 6Tr^^t,Lu,x)idy) F{duj) (E> dt - a.e. 

We fix now {t, w) G Q'j, and C, D G ^(M'^). 
We have by definition of T*'"^'" that 



^cninr],U^)^DiT''^^''ix))Ptidx), 



the latter because ^[r2'^](( j = 0. Notice that on the other hand, by definition of ^4*'", 

and T", for all (t, x) £ A*'" n fir we have that 

{y : (x,y) G 5npp(7r*'-'")} = = {T^{t,u;,x)}. 

This implies that C {(t,a;,x) G $7^ x M*^ : {y : (x,y) G snpp(7r*''^'")} is a singleton}. 
Now, let F*'"^ G ^(M"') be a measurable set with Pt(F*''^) = 1 and such that r*''^'''(x) = 
V(p*''^'''{x) is defined for all x G F*''^. Then, on F*'"^ n [j^T](t,a;) it must hold that 

r"(t,cj,x) = r*'^'"(x) = Vv9*'^'"(x), 

and we conclude that for all (t, tu) G 0^, 

ic(x)i^(r"(t,w,x))Pt((ix) 

□ 

We point out that Theorem 11.21 implies 

Corollary 5.3. r"(t,u;,x) = r*''^'"(x) holds ¥{dLj) O O Pt{dx)-almost surely . Conse- 
quently, r*''^'"(x) is measurable with respect to the completed a— field ofVred^ ^B{W^) with 
respect to ¥[du) ® dt® Pt{dx). 
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6 Pathwise convergence and rates for stochastic particle sys- 
tems to Landau process 

Proof of Proposition 14.31 

From the proof of Theorem 11.21 it is clear that integrals with respect to the measures 
l^k,n^Q^ Pt{dx) dt and lj^k,nPt{dx) dt are indistinguishable. By considering quadratic 
variations, the same is seen to hold for the stochastic integrals with respect to 1 j^u,n^^^Wp{dx, 
and lj^k,nWp{dx,dt). Write 

for the m— th coordinate of the process B^^'^ in ([7]), which is a real valued continuous local 
martingale with respect to J-^ (see [20]). Then, we have that 

=n%,™),(»',w) / / lAk,nr,A^',^r^^{s,uj,x)Ps{dx)ds 
Jo Jw 

= n'du^k,m),(i\k',m') / 'i-Ak,"nnrr(^^^^^)^s{dx)ds, 
Jo JR'i 

by step (a) in the proof of Theorem [L2j Now, for (s, uj) £ the points Xl{(jj), . . . , X'^{u;) 
are all different, and consequently we have that 

/ l^...nf7r(s,u;,x)P,((ix) =Ps{{x : r"(s, u;, x) = X,^(cu)}) 

=P,({x : T^'--"(x) = X^{u)}) 
=vr^'-'"({(x,y):y = X,^H}) 

1 

n 



Thus, we have l^B^^'^'^^B'^'''' '^''^')t = t5(^i^k,ni),{i',k',m')j and the result follows. 
We now are ready to prove Theorem ll.il 

Proof of Theorem 11.11, a) Let us fix n € N*, and define for i = 1, . . . n, 

rt n „t n 



□ 



or equivalently, in an indistinguishable way, 

xi' 



Jo J^'^t^i 

+ f [ YbiXl'^-X^ni^,,nis,y)Psidy)ds 
Jo JR"^ 
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By standard arguments and the fact that the sets A'''"' are disjoint (step (a) of the proof 
of Theorem 11.21) , we have 



E{\xr-xi\']< 




JM.'i 



\k=i ^ 



lj^k,n{s,y) Psidy)ds 



y) l^fe,n(s,y) ) I Ps{dy)ds 
(9) 



The first term in the right hand side of ([9]) is bounded by 



X^") lA.As,y)]\Ps{dy)ds 



1 2 



X^-)-a(X:-T-(s,y)) \ Ps{dy)ds 



JR'* 



+ ^/ / E\Y,[Hxi-T^{s,y))-a{Xi-y)Ylj,,,„{s,y))]Ps{dy)ds 



\k=l 



< 



^Ai',"is,y)Ps{dy) ds 



Jo Vfe=i ^ -^K' 
C j'^ j^^E(^{\T^{s,u:,y)-y\HA.As,y))^ Ps{dy)ds 



C / ^ Xr-X! ds 



+ C f I E(\T^{s,LO,y)-y\'')pAdy)ds 
Jo iR'* ^ ^ 

(\T^{s,u;,y)-y\^)Ps{dy)ds 



ds 




JR'' 
t 



= '^Cj E [jX'/' - Xl\ j ds + C j E{W^{u^,Ps))ds 

by exchangeabihty of ((X-*^'", X^), . . . , (X"'", X")). A similar bound is obtained for the 
second term in @. We deduce by Gronwah's lemma that 

E (|X,*'" - ) < Cexp(C'r) 1^ E{Wi {v^ , P,)) ds 
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By a little finer argument using a Burkholder-Davis-Gundy inequality, we can obtain as 
usual an estimate of the form 

E I sup \xi''' - Xl\A < Cexp{C'T) E{Wi{u^,Ps)) ds 
\te[Q,T] J Jo 

□ 

We recall a result proved in Rachev and Riischendorf [T3] giving L^-rates of convergence of 
empirical measures in the Wasserstein metric. 

Theorem 6.1. ( fU^ Theorem 10.2.1) Let ^ a probability on M'^ and let ,Y'^ , . . . ,¥"■ be 
independent identically distributed random variables with law /x. Let /x„ be the empirical 
measure of these variables. Then, if // has high enough finite absolute moments: c := 
/jjd \y\'^~^^ IJ-{dy) < oo, there is a constant C depending only on c and on the dimension d, 
such that 

E{wi{fXn,li)) <Cn^K 

Denote by W2 the Wasserstein distance between probability measures Q on the path space 

Ct ■■= C([0,r],R'^), such that /^^ supo<i<r \x{t)\'^Q{dx) < 00. 

From the previous result and Lemma l4.H it is simple to deduce the following 

Corollary 6.2. Let P be the pathwise law of the nonlinear process ([T|). Under the assump- 
tions of Theorem ] 1.1\ and moreover that J^a \y\'^~^^ Po{dy) < 00, we have that 

>V|(^au;(X^'"),P) < CT,dn^. 

□ 

The previous results are the first convergence rates obtained so far for stochastic particle 
systems of the "Landau type" ([1]), and they are not specific to the particular coefficients of 
the Landau equation ([2]). They justify the interest of the particle systems introduced in Q 
and are the first step in the construction and the numerical study of a simulation algorithm 
for {Pt)t- We notice that since we deal with space-time random fields, the dependence of 
the results on the dimension d is somewhat expectable, as opposite to the situation in the 
McKean-Vlasov model. The techniques we have introduced provide some insight about 
that dependence. 

Acknowledgements The authors are very grateful to Roberto Cominetti for helpful sug- 
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